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The spectral gap—the energy difference between the ground state and first excited state—is central to quantum many-body physics. 
Many challenging open problems, such as the Haldane conjecture, existence of gapped topological spin liquid phases, or the Yang-Mills 
gap conjecture, concern spectral gaps. These and other problems are particular cases of the general spectral gap problem: given a 
quantum many-body Hamiltonian, is it gapped or gapless? Here we prove that this is an undecidable problem. We construct families of 
quantum spin systems on a 2D lattice with translationally-invariant, nearest-neighbour interactions for which the spectral gap problem 
is undecidable. This result extends to undecidability of other low energy properties, such as existence of algebraically decaying ground- 
state correlations. The proof combines Hamiltonian complexity techniques with aperiodic tilings, to construct a Hamiltonian whose 
ground state encodes the evolution of a quantum phase-estimation algorithm followed by a universal Turing Machine. The spectral 
gap depends on the outcome of the corresponding Halting Problem. Our result implies that there exists no algorithm to determine 
whether an arbitrary model is gapped or gapless. It also implies that there exist models for which the presence or absence of a spectral 
gap is independent of the axioms of mathematics. 


The spectral gap is one of the most important physical prop¬ 
erties of a quantum many-body system, determining much of 
its low-energy physics. Gapped systems exhibit non-critical 
behaviour (e.g. massive excitations, short-range correlations). 
Whereas phase transitions occur when the spectral gap van¬ 
ishes and the system exhibits critical behaviour (e.g. massless 
excitations, long-range correlations). Many seminal results 
in condensed matter theory prove that specific systems are 
gapped or gapless. For example, Lieb, Schultz and Mattis’ 
proof that the Heisenberg chain is gapless for half-integer spin^ 
(later extended to higher dimensions by Hastings^), or Affleck 
et al.’s proof that the ID AKLT model is gapped.^ Similarly, 
many famous and long-standing open problems in theoreti¬ 
cal physics concern the presence or absence of a spectral gap. 
A paradigmatic example is the antiferromagnetic Heisenberg 
model in ID with integer spins. The “Haldane conjecture” 
that this model is gapped, first formulated in 1983,"^ has so far 
resisted rigorous proof despite strong supporting numerical 
evidence.^ The same question in the case of 2D non-bipartite 
lattices such as the Kagome lattice was posed by Anderson in 
1973.^ The latest numerical evidence' strongly indicates that 
these systems may be topological spin liquids. This problem 
has attracted significant attention recently^ as materials such 
as herbertsmithite^ have emerged whose interactions are well- 
approximated by the Heisenberg coupling. The presence of a 
spectral gap in these models remains one of the main unsolved 
questions concerning the long-sought topological spin liquid 
phase. In the related setting of quantum field theory, one of 
the most notorious open problems again concerns a spectral 
gap: the Yang-Mills mass gap problem.Proving a gap could 
provide a fuff explanation of the phenomenon of quark confine¬ 
ment. Whilst there is strong supporting evidence from lattice 
QCD numerics, the problem remains open. 

All of these problems are specific instances of the general 
spectral gap problem: Given a quantum many-body Hamilto¬ 
nian, is the system it describes gapped or gapless? Our main 


result is to prove that the spectral gap problem is undecidable in 
general. This says more than merely showing that the problem 
is computationally or mathematically hard. Though one may 
be able to solve the spectral gap problem in specific cases, our 
result implies that it is in general logically impossible to say 
whether a system is gapped or gapless. This has two meanings, 
and we prove both: 

1. The spectral gap problem is algorithmically undecidable: 
there cannot exist any algorithm which, given a description 
of the local interactions, determines whether the resulting 
model is gapped or gapless. This is the same sense in 
which the Halting Problem is undecidable.^^ 

2. The spectral gap problem is axiomatically independent: 
given any consistent recursive axiomatisation of mathe¬ 
matics, there exist particular quantum many-body Hamil¬ 
tonians for which the presence or absence of the spectral 
gap is not determined by the axioms of mathematics. This 
is the form of undecidability encountered in GodeTs in¬ 
completeness theorem. 

Precise statement of results 

It is important to be precise in what we mean by the spectral gap 
problem. To this end, we must first specify the systems we are 
considering. Since we are proving undecidability, the simpler 
the system, the stronger the result. We restrict ourselves to 
nearest-neighbour, translationally-invariant spin lattice models 
on a 2D square lattice of size L x L (which we later take to oo), 
with local Hilbert space dimension d. Any such Hamiltonian 
Hl is completely specified by at most three finite-dimensional 
Hermitian matrices describing the local interactions of the 
system: two d^ x d^ matrices /zj-ow and h^oi specifying the 
interactions along the rows and columns, and a J x J matrix hi 
specifying any on-site interaction. All matrix elements will be 
algebraic numbers, and we normalise the interaction strength 
max{||/zj-owll? Il^coill? 11^1 III — 1- 
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FIG. 1. Gapped and gapless systems, (a) A gapped system has 
a unique ground state Xq{H) and a constant lower-bound y on the 
spectral gap A(//) in the thermodynamic limit, (b) A gapless system 
has continuous spectrum above the ground state in the thermodynamic 
limit. 


(Hi). D is Hermitian, with matrix elements in {0, \,p}. 

For each positive integer n, define the local interactions of a 
translationally-invariant, nearest-neighbour Hamiltonian H(n) 
on a 2D square lattice as: 

h\{n) = a{n)IV 
hrow ~ 12 

hcoiin) = A+I3 C + , 

((pin) = is the rational number whose binary fraction 

expansion contains the binary digits of n after the decimal 
point, \(p\ denotes the number of digits in this expansion, a(n) 
is an algebraic number < fi computable from n, and 11 is a 
projector.) Then: 

(i). The local interaction strength is < 1 (i.e. ||/ti(^)||, ||/trowlL 

\\hcoiin)\\ < \). 



We must also be precise in what we mean by “gapped” and 
“gapless” (see Figure 1). Since quantum phase transitions occur 
in the thermodynamic limit of arbitrarily large system size, we 
are interested in the spectral gap 2s(Hl) = 2i(Hl) - Ao(Hl) 
as the system size L ^ oo (where Ao,Ai are the lowest and 
second-lowest eigenvalues). We take “gapped” to mean the 
system has a unique ground state and a constant lower bound 
on the spectral gap: > y > 0 for all sufficiently large L. 

We take “gapless” to mean the system has continuous spectrum 
above the ground state in the thermodynamic limit. 

Note that gapped is not the negation of gapless here; there 
are systems that fall into neither category. We adopt such 
strong definitions to deliberately exclude ambiguous cases, 
such as systems with degenerate ground states. A Hamiltonian 
that is gapped or gapless according to the above definitions 
is recognised as such throughout the literature. We show that 
the spectral gap problem is undecidable even with the promise 
that the Hamiltonian either has a unique ground state and a 
spectral gap of magnitude 1, or has continuous spectrum above 
the ground state. 

We prove this by showing that the Halting Problem for Tur¬ 
ing Machines can be encoded in the spectral gap problem, 
implying that the latter is at least as hard as the former. Recall 
that a Turing Machine is a simple, abstract model of compu¬ 
tation in which a head reads and writes symbols from some 
finite alphabet on an infinite tape and moves left or right, fol¬ 
lowing a finite set of rules. The Halting Problem asks: given 
an initial input written on the tape, does the Turing Machine 
halt on that input? Turing famously proved that this problem is 
undecidable,^^ and we relate it to the spectral gap problem in 
the following way: 

Theorem 1 We can construct explicitly a dimension d, df x df 
matrices A, B,C,D and a rational number fi > 0, which can be 
chosen as small as desired, so that: 

(i). A is Hermitian, with matrix elements in 7L + fiTL + 


(ii). If the universal Turing Machine halts on input n, then 
the Hamiltonian H(n) is gapped with y > 1. 

(in). If the universal Turing Machine does not halt on input n, 
then the Hamiltonian H(n) is gapless (i.e. has continuous 
spectrum). 

This implies that the spectral gap problem is algorithmically 
undecidable since the Halting Problem is. By a standard 
argument^'^ this also implies axiomatic independence. Both 
forms of undecidability extend to other low-temperature prop¬ 
erties, such as critical correlations in the ground state. In fact, 
our method allows us to prove undecidability of any physi¬ 
cal property that distinguishes a Hamiltonian from a gapped 
system with unique, product ground state. 

Hamiltonian construction 

We will first relate undecidability of the spectral gap to unde¬ 
cidability of another important physical quantity, the ground 
state energy density, which for a 2D lattice is given by 
Ep = limL^oo AoiHifilL^. We then show how to transform 
the Halting Problem into a question about ground state energy 
densities. 

Reducing the ground state energy density problem to the 
spectral gap problem requires two ingredients: 

(i). A translationally-invariant Hamiltonian F^m(^) on a 
2D square lattice with local interactions hu((p), whose 
ground state energy density is either strictly positive 
or tends to 0 from below in the thermodynamic limit, 
depending on the value of an external parameter (p. How¬ 
ever, determining which case holds should be undecid¬ 
able. Constructing such a Hamiltonian constitutes the 
main technical work of our result. (Note that these prop¬ 
erties are unaffected if we multiply hu((p) by an arbitrary 
fixed rational number fi, no matter how small.) 
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(ii). A gapless Hamiltonian with translationally-invariant 
local interactions hd and ground state energy 0.* There 
are many well-known examples of such Hamiltonians, 
e.g. the critical XY-model.^ 

Given Hamiltonians with these properties, we construct a new 
translationally-invariant Hamiltonian with local interactions 
h{ip) which is gapped or gapless, depending on the value of 
ip, in the following way. The local Hilbert space of h{(p) is 
the tensor product of those of hu and hd together with one 
additional energy level: 7T = |0) 0 0 TTj. We then take the 

interaction between nearest-neighbour sites i and j to be: 

= |O>(O|<'^0(1 -10>(0|)<^'^ 

( 1 ) 

Its spectrum can be seen to be (see supplementary information 
for details): 

spec H{ip) = {0} U I spec Hu{<p) + spec US, (2) 

with S > 1. Recalling that we chose Hd to be gapless, we see 
immediately from (2) that if the ground state energy density of 
Hu tends to zero from below (so that Aq(Hu) < 0) then H(cp) 
is gapless, but if Hu has strictly positive ground state energy 
density (so that X^iHu) diverges to +oo) then it has spectral gap 
> 1, as required (see Figure 2). 

Observe that this construction is rather general: by choosing 
different hd we obtain undecidability of any physical property 
that distinguishes a Hamiltonian from a gapped system with 
unique product ground state. 

Encoding computation in ground states 

To construct the Hamiltonian Hu{,<p), we encode the Halting 
Problem into the local interactions hu{<p) of the Hamiltonian. 
The latter concerns the dynamics of a classical system—a Tur¬ 
ing Machine. How can this be related to the ground state en¬ 
ergy density—a static property of a quantum system? The idea, 
which dates back to Feynman,^^ is to construct a Hamiltonian 
whose ground state encodes the entire history of the compu¬ 
tation in superposition. I.e. if the state of the computation at 
time t is represented by the state vector \il/t), the ground state is 
the so-called computational history state YJJq 10 In 
the following, we will sometimes refer to the Turing Machines 
encoded in the Hamiltonian “running” on some input. This 
should be understood to mean that the evolution produced by 
running the Turing Machine on that input appears in the ground 
state as the corresponding computational history state. 

If there are no other constraints, writing down such a Hamil¬ 
tonian is straightforward. However, constructing such a Hamil¬ 
tonian out of the local interactions of a many-body system is 
more involved. Feynman’s idea was substantially developed by 
Kitaev,^^ and after a long sequence of results^^"^^ culminated 
in Gottesman and Irani’s construction for ID spin chains with 
translationally-invariant, nearest-neighbour interactions.^^ 


* Recall from above that by “gapless” we mean continuous spectrum above the 
ground state, not merely a vanishing spectral gap. 
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FIG. 2. Ground state energy density to spectral gap. To relate ground 
state energy density and spectral gap, we need: (i) a Hamiltonian 
Hu((p) whose ground state energy density is either strictly positive 
or tends to 0 from below in the thermodynamic limit, but determin¬ 
ing which is undecidable, and (ii) a gapless Hamiltonian Hd with 
ground state energy 0. We combine Hu{<p) and Hd to form a new local 
interaction, h{<p), in such a way that h{<p) has (iii) an additional non¬ 
degenerate 0-energy eigenstate, and the continuous spectrum of Hd is 
shifted immediately above the ground state energy of H^. (a) If the 
ground state energy density of Hu{<p) tends to 0 from below, then its 
ground state energy in the thermodynamic limit must be < 0, and h{p) 
is gapless, (b) Whereas if the ground state energy density of Hu{<p) is 
strictly positive, its ground state energy in the thermodynamic limit 
must diverge to +oo, and h{<p) is gapped. 


For any Quantum Turing Machine^^ (QTM), Gottesman and 
Irani construct an interaction h between neighbouring particles, 
such that the ground state of the ID translationally-invariant 
Hamiltonian Hqi = Ki+i is of the form 10 Wt). 

where the “clock” part of the computational history state \t) ^ 
|l) 0 f counts time in unary, and represents the state 

of the QTM after t time-steps. Moreover, the ground state 
energy can be taken equal to 0. 

The translationally-invariant Hamiltonians we are consid¬ 
ering are completely specified by the finite number of matrix 
elements in the local interactions h^^^, h^ox and hi. To encode 
the Halting Problem in the Hamiltonian, we need a way of 
encoding any of the countably infinite possible inputs to the 
universal Turing Machine into these parameters. For these, we 
make use of quantum phase estimation.^^ 
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Quantum phase estimation 

Given a unitary U, the quantum phase estimation algorithm 
estimates an eigenvalue of (7 to a given number of bits 
of precision (which must be chosen in advance). It is well- 
known^^ that if the number of bits of precision in the quantum 
phase estimation algorithm is greater than or equal to the num¬ 
ber of digits in the binary fraction expansion of (p, then the 
quantum phase estimation algorithm, rather than estimating 
the phase approximately, will output all the digits of ip (written 
as a binary fraction) exactly. 

We use this property to construct a family of Quantum Tur¬ 
ing Machines indexed by ^ e N, with the following prop¬ 
erties: (i). The number of internal states and tape symbols of 
Pn are independent of n. (ii). Given a number N = 2^ -\>n 
as input (written in binary), Pn writes the binary expansion 
of n on its tape and then halts deterministically. (The reason 
for having input N of this form will become clear later). To 
construct we construct a QTM which uses the input N 
to determine how many digits of precision to use, then runs 
the quantum phase estimation algorithm on the single-qubit 
gate (/ = )• The phase ^ in 1/ is determined by the 

transition rules of the QTM.^^ Choosing ip to be the rational 
number whose binary fraction expansion contains the digits of 
n (expressed in binary) achieves the desired behaviour for 
By “dovetailing” Pn with a universal Turing Machine UTM 
(i.e. run Pn first, then run UTM), the UTM runs on the input 
specified by ip. 

The quantum computation carried out by Pn followed by 
UTM can be encoded in the Hamiltonian using the history 
state construction described above. The phase being 
estimated then becomes one of the matrix elements of the 
Hamiltonian. The same happens with the appearing in 
the inverse Quantum Fourier Transform—the key ingredient 
of the quantum phase estimation algorithm. 

Finally, we must ensure that the |0^o) component of the his¬ 
tory state is correctly initialised to input of the form = 2^ - 1 
(written in binary) required by But = 2^ - 1 in binary is 
simply a string of A^ H’s, and it is easy to ensure that |</ro) is the 
state |1)®^ using translationally-invariant local interactions. 

If we add an on-site interaction h\ - |T)(T| to the history 
state Hamiltonian constructed above, giving an additional en¬ 
ergy to the state |T) representing the halting state, then its 
ground state will pick up additional energy if and only if the 
UTM halts. However, the ground state energy still converges 
to 0 with system size in both cases (see supplementary infor¬ 
mation). The energy density therefore tends to zero in the 
thermodynamic limit, whether or not the UTM halts. 

To remedy this, and amplify the difference between the halt¬ 
ing and non-halting cases, we use the second spatial dimension 
and exploit Wang tilings. 

Quasi-periodic tilings 

A Wang tile^^ is a square with markings along each edge. A 
tiling is then an arrangement of such tiles covering the whole 
plane, so that the markings on adjacent edges match. A tiling 
can easily be encoded in a ground state of a classical Hamil¬ 
tonian on a 2D square lattice: by representing tile types by 
an orthogonal basis (IT)} for the local Hilbert space 77^ and 


choosing local interaction terms \Ti){Ti\ 0 \Tj){Tj\ to give an 
energy penalty to all adjacent non-matching pairs of tiles Tt, Tj, 
a tiling of the plane is equivalent to a ground state with energy 
zero. 

We prove and exploit very particular properties of an aperi¬ 
odic tiling due to Robinson,^^ and combine this with the history 
state Hamiltonian. Although the pattern of tiles in the Robin¬ 
son tiling extends infinitely in all directions, it never repeats. 
More precisely, it contains periodically repeating sub-patterns 
forming squares with sizes given by every power of 4 (see 
Figure 3). This will allow us to encode in the ground state 
many copies of the UTM running on the same input ip, with 
tapes of all possible finite lengths and for every possible finite 
run-time (see Figure 3). 

This is achieved by sandwiching the ID quantum history- 
state Hamiltonian hq “on top of” the Robinson tiling Hamil¬ 
tonian he to form two “layers”, so that the local Hilbert space 
at each site is 7T = 7-4 0 (7-4 © "Hq) (where 7-4 = |0) is an 
additional energy level). One can then construct a Hamiltonian 
(see supplementary information) whose ground state is of the 
form |r)c 0 \^)eq^ where \T)e is a. product state representing 
a classical configuration of the tiling layer and contains 
- in a tensor product structure - computational history states 
along one edge (called a “segment”) of all of the squares ap¬ 
pearing in the configuration given by T. These computational 
history states are essentially the only constituents of which 
contribute to the energy. The Hamiltonian also has an on-site 
interaction hi = |T)(T| that gives an additional energy to the 
state |T) representing the halting state of the Turing Machine. 
Hence the ground state will pick up additional energy from all 
encoded Turing Machines that halt. This energy still decreases 
with the effective size of the system which, however, is now 
the size of the corresponding segment in the Robinson tiling 
(see Figure 3), not the overall system size. 

What can we say about the ground state energy? If the UTM 
does not halt on input n, then \T)^ is a valid tiling and for 
all segments which are larger than \n\ the ground state energy 
contribution is 0. The contribution for each segment smaller 
than 1^1 is some algebraic computable number. If a(n) is the 
constant that adds up the contributions of all segments smaller 
then 1^1, the addition of the energy shift hi = -a(n)l to the 
Hamiltonian makes the ground state energy density negative 
(and tending to 0 from below) in the non-halting case (see 
supplementary information). 

In the halting case, one of two things may happen. If \T}e 
is a valid tiling, the number of squares large enough for the 
encoded Turing Machine to halt grows quadratically with sys¬ 
tem size, and each of them contributes a small but non-zero 
energy. Since such a state also picks up the energy contribu¬ 
tion from segments of size smaller than \n\, even after adding 
hi = -a(n)l the energy diverges with lattice size. Hence the 
ground state energy density is strictly positive in the halting 
case, as desired. 

Alternatively, one could try to reduce the energy by intro¬ 
ducing defects in the tiling, that effectively “break” some of 
the Turing Machines so that they do not halt. However, we 
prove that the Robinson tiling is robust against such defects: 
a tile mismatch can only affect the pattern of squares in a fi- 
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FIG. 3. Complete Hamiltonian construction. The Robinson tiles 
enforce a recursive pattern of interlocking squares, of sizes given by 
all integer powers of 4 (middle figure). As with any Wang tiling, we 
can readily represent this as a classical Hamiltonian whose ground 
state has the same quasi-periodic structure. Since the set of tiles is 
fixed, the local dimension of this Hamiltonian is constant. By adding 
a “quantum layer” on top of this tiling Hamiltonian and choosing a 
suitable translationally-invariant coupling between the layers, we can 
effectively place copies of QTM encoded in a ID history state Hamil¬ 
tonian (top figure) along one edge of all of the squares. The ground 
state of this Hamiltonian consists of the Robinson tiling configuration 
in the tiling layer, with computational history states in the quantum 
layer along one edge of all of each square in the tiling (bottom figure). 
Each of these encodes the evolution of the same quantum phase esti¬ 
mation algorithm and universal Turing Machine. The effective tape 
length available each QTM is determined by the size of the square it 
“runs” on. 


nite region around the defect, and each defect contributes 0(1) 
energy. We can choose the parameters (see supplementary in¬ 
formation) to guarantee that introducing defects is energetically 
unfavourable. This completes the argument for our main result. 
Additional technical details can be found in the supplementary 
information (online). 

Discussion 

Let us take a step back and discuss both the implications and 
the limitations of this result. The result clearly concerns math¬ 


ematical models of quantum many-body systems, but also the 
behaviour of and methods for treating the thermodynamic limit. 
Moreover, the result can also be seen as an indication of new 
physical phenomena. 

An immediate consequence of our result is that there cannot 
exist an algorithm or a computable criterion that solves the 
spectral gap problem in general. Whilst algorithmic undecid¬ 
ability always concerns infinite families of systems, the ax¬ 
iomatic interpretation of the result also allows us to talk about 
individual systems: there are particular Hamiltonians within 
these families for which one can neither prove nor disprove 
the presence of a gap, or of any other undecidable property. 
Unfortunately, our methods cannot pinpoint these cases, let 
alone prove that one of the big open problems mentioned in 
the introduction is of this kind. 

A further consequence concerns the behaviour of the ther¬ 
modynamic limit. In practice, we usually probe the idealised 
infinite thermodynamic limit by studying how the system be¬ 
haves as we take finite systems of increasing size. One often 
assumes that the systems, though finite, are so large that we 
already see the asymptotic behaviour. In numerical simulations 
of condensed matter systems, one typically simulates finite 
systems of increasing size and extrapolates the asymptotic 
behaviour from finite size scaling. Similarly, lattice QCD 
calculations simulate finite lattice spacings, and extrapolate the 
results to the continuum. Renormalisation group techniques 
accomplish something similar mathematically.^^ However, the 
undecidable quantum many-body models constructed in this 
work exhibit behaviour that defeats such approaches. As the 
system size increases, the Hamiltonian will initially look like 
a gapless system, with the low-energy spectrum appearing to 
converge to a continuum. But at some threshold lattice size, 
a spectral gap of magnitude 1 will suddenly appear (or, vice 
versa, a gap will suddenly close^^). Not only can the lattice 
size at which the system switches from gapless to gapped be 
arbitrarily large, the threshold at which this transition occurs 
is uncomputable. The analogous implication also holds for 
all other undecidable low-temperature properties. Thus any 
method of extrapolating the asymptotic behaviour from finite 
system sizes must fail in general. 

This leads us directly to new physical phenomena. First, it 
hints at a new type of “phase transition”, not driven by tempera¬ 
ture or extrinsic local parameters, but by the size of the system. 
Some of the models constructed in the proof exhibit a drastic 
and abrupt change of properties when their size is increased 
beyond a certain scale. The scale at which this happens can be 
very large and is, in fact, not generally computable from the 
local description of the system. Second, our results show that 
certain quantum many-body models exhibit a radical form of 
instability. An arbitrarily small change in the parameters can 
make the system cross an arbitrary number of gapped/gapless 
transitions. In a sense, this phenomenon is the source of the 
undecidability in our models. 

We finish this discussion with a closer look at some of the 
limitations of the result. First, our results concern 2D (or 
higher-dimensional) systems. Although the majority of our 
construction is already ID, we do not currently know whether 
the entire result holds in ID as well. Second, although a 
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theoretical model of a quantum many-body system is always 
an idealisation of the real physics, the models we construct 
in the proof are highly artificial. Whether the results can be 
extended to more natural models is open. A related point is that 
we prove our results for Hamiltonians with a very particular 
form. We do not know how stable the results are to small 
deviations from this. This is a general issue with most many- 
body models; stability in this sense is not understood even 
for much simpler models such as the Ising model. Recent 
stability proofs only apply to certain types of frustration-free 


Hamiltonian. Indeed, our results restrict the extent to which 
such stability results can be generalised. Similarly, we do 
not know whether the results hold for systems with small 
local Hilbert space dimension. Whilst the dimension d in 
Theorem 1 is fixed and finite, providing an estimate for it 
would be cumbersome and certainly involve large exponentials. 
However, the steps in the proof described above are not tailored 
to minimising this dimension. Whether there is a non-trivial 
bound on the local Hilbert space dimension below which the 
spectral gap problem becomes decidable is an intriguing open 
question. 
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